Chaos in Quantum Dots: Dynamical Modulation of Coulomb Blockade Peak Heights 
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We develop a semiclassical theory of Coulomb blockade 
peak heights in quantum dots and show that the dynamics in 
the dot leads to a large modulation of the peak height. The 
corrections to the standard statistical theory of peak height 
distributions, power spectra, and correlation functions are 
non-universal and can be expressed in terms of the classical 
periodic orbits of the dot that are well coupled to the leads. 
The resulting correlation function oscillates as a function of 
peak number in a way defined by such orbits; in addition, the 
correlation of adjacent conductance peaks is enhanced. Both 
of these effects are in agreement with recent experiments. 

PACS numbers: 73.23.Hk, 05.45.Mt, 73.20.Dx, 73.40. Gk 

The electrostatic energy of an additional electron on 
a quantum dot- a mesoscopic island of confined charge 
with quantized states- blocks the flow of current through 
the dot, an effect known as the Coulomb blockade |Q. 
Current can flow only if two different charge states of 
the quantum dot are tuned to have the same energy; 
this produces a peak in the the conductance of the dot 
whose magnitude is directly related to the magnitude of 
the wavefunction near the contacts to the dot. Since 
dots are generally irregular in shape, the dynamics of the 
electrons is chaotic, and the characteristics of Coulomb 
blockade peaks reflect those of wavefunctions in chaotic 
systems Previously, a statistical theory for the 

peaks was derived by assuming these wavefunctions 
to be completely random and uncorrelated with each 
other. The experimental data f° r the distribution of 
the Coulomb blockade peak heights were found to be in 
excellent agreement with the predictions of the statistical 
theory, thus giving a solid foundation to the conjecture 
of effective "randomness" of the quantum dot wavefunc- 
tions. 

It therefore came as a surprise when several recent ex- 
periments [^|-^) demonstrated large correlations between 
the heights of adjacent peaks. The effect of nonzero tem- 
perature (when several resonances contribute to the same 
peak) was found to be insufficient to account for these 
correlations Q . To explain the correlations, the enhance- 
ment due to spin-paired levels , due to a decrease of 
the effective level spacing found in density functional cal- 
culations Q , and due to of level anticrossings in inter- 
acting many-particle systems pd| ] were proposed. How- 
ever, we show here that peak height correlations arise 
already within an effective single-particle picture of the 
electrons in the quantum dot. The specific internal dy- 
namics of the dot, even though it is chaotic, modulates 



the peaks: because all systems have short-time features, 
chaos is not equivalent to randomness. The predicted 
dynamical modulation is exactly of the type needed to 
explain the experiments |^-^| ■ 

To study the non-universal effects of the dynamics of a 
particular dot, we derive a relation between the quantum 
conductance peak height and the classical periodic orbits 
in the dot. The main effect is that as a system parameter 
varies- the magnetic field or the number of electrons in 
the dot in response to a gate voltage, for instance- the in- 
terference around each periodic orbit oscillates between 
being destructive and constructive. When the interfer- 
ence is constructive for those periodic orbits which come 
close to the leads used to contact the dot, the wavefunc- 
tion is enhanced near the leads, the dot-lead coupling 
is stronger, and so the conductance is larger. Likewise, 
destructive interference produces a smaller conductance. 
The resulting modulation can be substantial, as shown in 
Fig. 1. Similar short-time dynamical effects have been 
noted in other contexts such as atomic and molecular 
spectra |T^-Q, eigenfunction scarring flf4] , |l~5| | , magne- 
totransport in antidot lattices Jl6| ], and tunneling into 
quantum wells Jl7[^19[. Such modulation is completely 
omitted in theories in which the wavefunction is assumed 
to change randomly as the system changes ^,|) . 
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FIG. 1. The peak conductance from tunneling through 
subsequent energy levels for the stadium billiard shown in 
the inset. Each peak is placed at the wavevector k corre- 
sponding to its level; 7? is the radius of the half-circle parts of 
the stadium. A Gaussian lead wavefunction appropriate for 
tunneling from a single transverse mode is used with width 
fea e ff = 15. 

Our starting point is the connection between the peak 
height and the width of the level in the quantum dot. We 
consider a dot close to two leads (Fig. I inset) so that the 
width, r, of a level comes from tunneling of the electron 
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to either lead. When the mean separation of levels is 
larger than the temperature T which itself is much larger 
than the mean width, the electrons pass through a single 
quantized level in the dot. The conductance peak height 
in this regime is p(J] 

r - e2?rr m 

pcak " h AkT { } 

Here for simplicity we consider symmetrically placed 
leads- the total width is equally split between tunnel- 
ing to the right and left leads- spinless particles, and 
temperature much smaller than the level spacing. 

The width of the level is related by Fermi's Golden 
Rule to the square of the matrix element for tunneling 
between the lead and the dot, M l ~* d . A convenient ex- 
pression for the matrix element in terms of the lead and 
dot wavefunctions, 4^ and respectively, was derived 
by Bardccn [El| and can be expressed as 



M ^ d = n_ r jsvttfvMr) 

m* Js 



(2) 



where the surface S is the edge of the quantum dot. T, 
then, depends on the square of the normal derivative of 
the dot wavefunction at the edge weighted by the lead 
wavefunction. Writing the product of the two ^d's in 
r as a Green function G(r, r') and using the standard 
semiclassical expression for the latter p2|, we express T 
as a sum over the classical trajectories which start at r 
on the edge of the dot near the lead and end at r' which 
is also on the edge near the lead. 

Tunneling from the lead to the dot is dominated by 
the lowest transverse energy subband in the constriction 
between the lead and the dot S. Therefore, for the cal- 
culation of the tunneling matrix clement the transverse 
potential in the tunneling region can be taken quadratic: 
Ue ~ k (y — yg) . In this case the transverse depen- 
dence of the lead wavefunction is simply a harmonic 
oscillator wavefunction, so that at the edge of the dot 
&t ~ Cfexp[— {y — ye) 2 /2a 2 s ], where y represents the 
direction tangential to the boundary of the dot, ye is 
the center of the lead and constriction, and the effective 
width is deft = Vh/ \j2nm* . While the exact form of 
the lead wavefunction is not crucial, the ^-dependence 
of the width is important for the semiclassical argument 
which follows; note that a c s ~ Vh does not depend on a 
particular transverse potential. 

Using this information about 4 e in the expression for 
M , we see that the lead wavefunction restricts the 
integration to a semiclassically narrow region of width 
a c g ~ Vh. This allows one to express the contribution 
of the open trajectories entering the Green function in 
terms of an expansion near their closed neighbors. In 
the resulting expression for T, the contribution of each of 
these closed orbits is suppressed by a factor exponentially 
small in Ap 2 , where Ap y is the change of transverse mo- 
mentum after one traversal. This suppression is the effect 



of the mismatch of the closed orbit (momentum) with the 
distribution of transverse momentum at the lead, which 
is centered at zero with width Spe ~ h/a c g ~ Vh for the 
lowest subband. Therefore, only closed orbits with semi- 
classically small momentum change Ap contribute to the 
width. This in turn implies that the closed orbit is lo- 
cated semiclassically close (within a distance ~ Vh) to a 
periodic orbit for which Ap = 0. Thus, one can express 
the tunneling width in terms of the properties of these 
periodic orbits, obtaining E2] 



r = r+ £ a, 



Ai:p.o. 

where the monotonic part is 



(3) 
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the amplitude is 



A, = 4V2 ^ [Tr 2 [M„] (1 + o-l) (1 + al )] ^ 
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v m 2 i + mi2, 
cos6» (y^ - ye) /a cff + sin6» p%a cS /h 
cos 6 p%a eS /h - sin (y^ - ye) /a cS , 



and, finally, the result for the slowly varying phase (j)^ 
will be given elsewhere. Here /„ is the Bessel function of 
complex argument, p M is the electron momentum for the 
periodic orbit ji at the bounce point, y^ is the bounce 
point coordinate, is the action of the periodic orbit, 
and Mfj, = (jnfj ) is the corresponding monodromy matrix 
p2| . The semiclassical approach used here is similar to 
the calculation of the tunneling current in a resonant 
tunneling diode in a magnetic field developed in Ref. fTij] 
and @. 

The equation above is the main result of this paper: it 
expresses the modulation of the heights of the Coulomb 
blockade peaks by the classical periodic orbits. Note that 
the result (||) is valid not only for chaotic but also for both 
integrable and mixed systems (for an integrable system 
or for the contributions of the remaining unbroken tori 
of a mixed system Tr [M] = 2). 

In order to assess the validity of the semiclassical ex- 
pression (|3|), we compare it to numerical calculations for 
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two simple billiards, one integrable- the circle- and one 
chaotic- the stadium. The stadium billiard (Fig. 1 inset) 
is one of the canonical examples of a completely chaotic 
system ]l2[] . From the dot wavefunctions (numerically 
obtained using the method of Rcf. [^3|), T can be calcu- 
lated from Eq. (Q) using the ^£ given above. To observe 
the variation in peak height, we vary the energy, or equiv- 
alently the wavevector k—p/h, which changes the number 
of electrons on the dot as more levels are filled. 
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FIG. 2. Length spectrum of the oscillations in G(k) ob- 
tained from the Fourier power, numerical (blue) and semiclas- 
sical (red) results compared. The power is normalized to the 
mean conductance and then this mean is removed for clar- 
ity, (a) Circular dot with narrow leads, ka e g ~ 1.2. where 
a e a is the width of the lead wavefunction. (b) Circular dot 
with wider leads, ka c s ~ 12. (c) Stadium dot using data in 
Fig. 1; dependence of amplitude at L/R= 4 on ka c g in in- 
set. The width of the peaks reflects the length of G(k) used. 
More data was available for the circle because it is integrable; 
conservation of angular momentum allows a simple represen- 
tation of the wavefunctions in terms of Bessel functions. In 
(a) the peak at L/R—4 is the diameter, that at 8 is its rep- 
etition, those accumulating to 2ir are the whispering gallery 
trajectories, and the largest length peak is the star orbit. The 
magnitude of the oscillatory part compared to the mean de- 
pends on the strength of the coupling to the periodic orbit 
and so on a e ff as well. For the stadium, (c), the principal 
peak corresponds to the horizontal orbit, which appears at 4 
because we use only the wavefunctions symmetric about the 
vertical symmetry axis (equivalent to using a half-stadium). 
Note the excellent agreement between the semiclassical theory 
and the numerical results in all cases. 
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FIG. 3. Conductance statistics: (a) peak-to-peak corre- 
lation function and (b) probability distribution of G for the 
stadium data in Fig. 1. The numerical correlation function 
(circles with typical error bars)- the average of all pairs of 
peaks m peaks apart- is in good agreement with the semiclas- 
sical theory (red). The agreement for small m is surprising 
since this regime is not semiclassical, but shows how dynamics 
can give rise to correlations even between nearest-neighbors. 
The numerical probability distribution (histogram) is for the 
entire range of data in Fig. 1 and is compared to both the 
semiclassical theory (red) and the standard statistical theory 
based on random wavefunctions (blue) . The two theories pre- 
dict nearly the same result for this quantity, and both are 
consistent with the numerics. 

In Fig. 1 we show an example trace of G for the sta- 
dium. The calculation clearly demonstrates both strong 
peak-to-peak fluctuations and an oscillatory modulation 
of the heights (3 periods are observed). While the former 
comes from the quasi-random fluctuations in the wave- 
functions near the leads, the large oscillatory modulation 
is caused by interference along the horizontal orbit which 
connects both leads. 

Since the main theoretical result concerns the periodic 
modulation of the peak heights, it is natural to consider 
the Fourier power spectrum of G pea k(k). In Fig. 2 we 
present a comparison of the numerical and semiclassi- 
cal power spectra, calculated for both integrable (circu- 
lar) and chaotic (stadium) dots. The data show that for 
both the circle and the stadium the power spectrum has 
sharp peaks corresponding to periodic orbits. More peaks 
appear for narrow leads [Fig. 2(a)] because the lack of 
momentum constraint in this regime allows coupling to 
more periodic orbits. The excellent agreement between 
the semiclassical expression and the numerical result in 
all cases is a striking demonstration of the validity of our 
theory. 

Further characterization of the peak fluctuations is 
shown in Fig. 3. The peak-to-peak correlation function 
is 



Cc 



[5G, SG] 



= (6G{E n+m )6G (E n )) n 
{SG{E n )X 



(4) 



where SG (E m ) = G (E m ) — (G (E n )) n is a natural mea- 
sure of the statistics of nearby peaks. A semiclassical 
expression for this quantity can be derived by assuming 
that the distribution of individual peak heights is locally 
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Porter-Thomas [ 24| , with the mean given by the semiclas- 
sical envelope (p). Indeed, as was first shown by Kaplan 
and Heller JL5], this is generally true for wavefunction 
fluctuations in chaotic systems. We obtain 



Corr„ 



<5m,0 + (1 - S m fi) X 



4r 2 + 3E u ^ 



(5) 



In Fig. 3(a) we compare the semiclassical correlation 
function with numerical data for the stadium dot. The 
oscillatory behavior for large separations reflects the peak 
in the corresponding power spectrum in Fig. 2 and is in 
agreement with the semiclassical result. The positive cor- 
relation for nearest neighbors is also in agreement with 
the semiclassical theory, demonstrating the influence 
of dynamics even in this apparently non-semiclassical 
regime. 

When T^> A, the major source of correlations between 
neighboring peaks is the joint contribution of several res- 
onances to the same conductance peak || . In this regime 
the "nearest- neighbor" correlator is Corr 7n= i~l, and the 
dynamical effect accounts for only a small correction to 
the correlation function. However, for low temperature 
T< A, the correlations due to temperature are exponen- 
tially suppressed. In this regime, the correlations induced 
by dynamical modulation dominate, and they account for 
the experimentally observed enhancement of correlations 
at low temperatures [|| . 

The probability distribution of G pca k over a large 
energy range is the main quantity considered in the 
previous statistical theories fill. They predict no 
peak-to-peak correlation or periodic modulation of the 
heights, and a Porter-Thomas distribution: -P(G pca k) = 
y/ 4/7rG pC ak exp(-Gp 0a k). Considering an energy range 
larger than any period in Eq. (^|), we find, in contrast, 
that the distribution should be locally Porter-Thomas 
but with the mean modulated by the periodic compo- 
nents, as in Ref. [ fl5| . Curiously, the resulting distri- 
bution is not very different from Porter-Thomas: Fig. 
3(b) shows that the two theories predict nearly the same 
result, and both are consistent with numerical calcula- 
tion. This explains why no dynamical effect was observed 
in the experimental peak-height probability distribution 

SB- 

In contrast, the periodic modulation of the peak 
heights has been observed in several recent experiments 
The clearest observation is in Ref. ||: the data in 
their Fig. 1 show modulated peak heights as a function 
of the number of electrons in the dot. In their trace of 90 
peaks, approximately six oscillations are visible, yielding 
a period of ~ 15 peaks. In our treatment, this period is 
the ratio of the period of fundamental oscillation in Eq. 
(||) to the level spacing A. The fundamental period is 

given by ^^§f^j = ?&A> where r M is the period of the 
relevant orbit. To determine t m , we use the billiard ap- 



proximation: = vpT^, where is the length of the 
shortest orbit and vf is the Fermi velocity. We use the 
micrograph of the dot to estimate for the T^-shaped 
orbit connecting the two leads, and calculate vp from the 
experimental density p5[ . Using the appropriate spin- 
resolved level spacing A= 10 [j,eV (which is half of the 
spin- full value from excitation measurements in Ref. 0), 
we find 7i/t^A«12. Because the billiard approximation 
underestimates the period in a soft wall potential, this is 
a lower bound for the modulation period, and therefore 
our theory is in good agreement with the experiment. 

Similarly, we make estimates which are consistent with 
the other two experiments showing variation as a func- 
tion of number of electrons [§J^]. A similar approach 
to the peak modulation as a function of magnetic field 
is in agreement with the experimental results [[||7j as 
well. This agreement with experiment is perhaps sur- 
prising, since the adding of electrons changes the effec- 
tive potential defining the dot; however, experiments on 
" magnetofingerprints" of the peaks 1 26 1 suggest that this 
change is small while to affect the dynamical modulation 
one must substantially change the action of the shortest 
periodic orbit. 

We close with two further experiments suggested by 
our results. First, if the tuning parameter used to change 
the number of electrons, such as a gate voltage, does not 
change the action of the dominant periodic orbit, then 
no modulation connected to that orbit should be seen. 
In particular, gates which affect different parts of the dot 
may produce different oscillatory behavior. Second, sev- 
eral samples made in a robust geometry- a circle with 
directly opposite leads, for example- should show the 
same modulation. Any deviations from the same behav- 
ior would be a sensitive indication of the material quality. 
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